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Abstract 

The hydrodynamic expansion of quark-gluon plasmas with spherical and longitu- 
dinally boost-invariant geometries is studied as a function of the initial energy density. 
The sensitivity of the collective flow pattern to uncertainties in the nuclear matter 
equation of state is explored. We concentrate on the effect of a possible finite width, 
AT ~ 0.1 T c , of the transition region between quark-gluon plasma and hadronic phase. 
Although slow deflagration solutions that act to stall the expansion do not exist for 
AT > 0.08 T c , we find, nevertheless, that the equation of state remains sufficiently soft 
in the transition region to delay the propagation of ordinary rarefaction waves for a 
considerable time. We compute the dependence of the pion-interferometry correlation 
function on AT, since this is the most promising observable for time-delayed expansion. 
The signature of time delay, proposed by Pratt and Bertsch, is an enhancement of the 
ratio of the inverse width of the pion correlation function in out-direction to that in 
side-direction. One of our main results is that this generic signature of quark-gluon 
plasma formation is rather robust to the uncertainties in the width of the transition 
region. Furthermore, for longitudinal boost-invariant geometries, the signal is likely to 
be maximized around RHIC energies, ~ 200 AGeV. 



' This work was supported by the Director, Office of Energy Research, Division of Nuclear Physics 
of the Office of High Energy and Nuclear Physics of the U.S. Department of Energy under Contract No. 
DE-FG-02-93ER-40764. 
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1 Introduction 



Lattice calculations of the thermodynamical functions of quantum chromodynamics (QCD) 
indicate ]l] that, (at zero net baryon density) in the vicinity of a critical temperature T c ~ 
160 MeV, strongly interacting matter undergoes a rapid transition from a (chirally broken, 
confined) hadronic phase to a (chirally symmetric, deconfined) quark-gluon plasma (QGP). 
The width of that transition region is presently only known to be in the range < AT < 
0.1 T c ~ 16 MeV. Therefore, one cannot yet conclude whether the transition is a first order 
phase transition (AT = 0), or merely a rapid increase of the entropy density associated with 
the change from dn hadronic to cIq quark and gluon degrees of freedom. 

One of the primary goals of relativistic heavy-ion physics is the creation and experimental 
observation of the predicted QGP phase of matter. Many signatures have been proposed such 
as electromagnetic radiation of thermal dileptons and photons 0, J/^f -suppression ||, jet 
quenching [f|], strangelet formation ||, or disordered chiral condensates (DCC's) 0. These 
signatures, however, do not depend directly on the actual form of the nuclear matter equation 
of state. Thermal electromagnetic radiation is, for instance, generic to any hot system, 
independent from its degrees of freedom (as long as they have electromagnetic charge). For 
example, it was shown @ that a hot hadron gas shines as brightly as a QGP. Similarly, 
jet quenching and J/\l/-suppression[] are generic consequences of final state interactions in 
any form of dense matter ||. Finally, strangelet or DCC formation require very specific 
assumptions about the dynamical evolution of the system. 

It is therefore of interest to study signals that are more directly related to the QCD 
equation of state. Signals of this type emerge from the influence of the equation of state 
on the collective dynamical evolution of the system. Such signals must be calculated within 
the framework of relativistic hydrodynamics ||10|| , since that is the only dynamical model 
which provides a direct link between collective observables and the equation of state. Of 
course, the use of such an approach requires a strong dynamical assumption, namely that 
the equilibration rates are much larger than typical gradients of thermodynamic quantities 
in the system. It now appears that radiative gluon energy loss in a QGP may be sufficiently 
large JO]] to support local equilibration on time scales less than 1 fm/c. In the following, 
we therefore neglect dissipative effects and assume the validity of ideal hydrodynamics to 
compute the collective evolution of the system. 



It was shown in [12, [L3L 14| that the transition to the QGP softens the equation of state 



in the transition region, and thus reduces the tendency of matter to expand on account of its 
internal pressure. This, in turn, delays the expansion and considerably prolongs the lifetime 



of the system. In ||15| , |16] it was moreover shown that this effect leads to a reduction of the 
transverse directed flow in semi-peripheral collisions that can be readily tested experimentally 
at fixed target energies [|l7j. 

In [13 the flow structure for the one-dimensional expansion of a slab of matter ( "Landau 
expansion model" |18[) was studied as a function of its initial energy density e and of the 



width AT of the transition region in the equation of state. In particular, it was shown how, 



Renewed interest in that signature, has, however, emerged with the advent of recent results from the 
NA50 experiment at the CERN SPS, showing a rapid variation of J/vE'-suppression as a function of transverse 
energy in Pb+Pb-collisions at 158 AGeV §. 



2 



in addition to the softening of the equation of state, for AT < AT* ~ 0.07676 T c the type of 
the hydrodynamical expansion solution changes qualitatively from a simple rarefaction wave 
to a rarefaction discontinuity. Such deflagration fronts |19| have small velocities, Vd <C c, 



and thus the conversion of QGP matter into hadronic matter is considerably delayed. In [II 
we also showed that for larger transition widths (AT > AT*) the lifetime of matter with 
temperature around T c is considerably reduced but cooler regions with T ~ 0.7 T c remain 
long-lived. 



In the present work we extend the investigations of [TJj] to more realistic, 3-dimensio- 
nal geometries. We consider systems with spherical symmetry ("fireball" geometry), and 
with transverse cylindrical symmetry and boost-invariant initial conditions in longitudinal 
direction ("Bjorken cylinder" geometry |2(|). For these symmetric geometries it is possible 
to use a modified one-dimensional algorithm as e.g. presented in fl3|]. Note that in contrast 



to the simple one-dimensional Landau geometry, single inclusive particle spectra and also 
two-particle correlation functions can be calculated in a realistic way from the hydrodynamic 
solutions for the above geometries. 

In this paper we show that the influence of the width AT of the transition on the 
dynamical evolution of the system in these 3-dimensional geometries is very similar to that 
found in [|TJ]] for purely one-dimensional expansion. In particular, the softening of the 
equation of state associated with a rapid cross-over region leads to a delay in the expansion 
and a prolonged lifetime of matter with temperature below T c . 

To explore observable consequences of the prolonged lifetime of the system, we then cal- 
culate the inverse width of the two-particle correlation functions in side- and out-direction, 
-R S ide and -R ou t- We demonstrate that the initial energy density dependence of the system's 
lifetime is mapped closely by the excitation function of the ratio -R ou t/-Rside- The phase tran- 
sition leads to a time-delayed expansion that is reflected in an enhancement of the above 
ratio relative to its value obtained in the ideal gas case without a transition. This signature 



of the QGP transition was first proposed by Pratt |2Tfl who performed similar calculations 
for the fireball geometry^ and, independently, by Bertsch et al. who employed a kinetic 
model for the hadronization in the Bjorken cylinder geometry. These previous investigations 
assumed the existence of a sharp first order transition, AT = 0. We show here that this 
signature is a generic feature of a rapid cross-over transition and holds for AT = 0.1 T c as 
well. 

The importance of the time delay as a signature of QGP formation is that it is among 
the few signatures that probe directly the equation of state of strongly interacting matter. It 



tests whether there exists a "soft region" in the equation of state of ultra-dense matter [F2 



which is most naturally characterized [14]] by a reduction of the local velocity of sound c$ as a 



function of the local energy density. As discussed in detail in the next section, a very general 
feature of a transition (that happens within a sufficiently narrow temperature interval AT) 
is the existence of such a region. For matter passing through that region of energy densities 
during the expansion phase, the flow will temporarily slow down or even possibly stall under 
suitable conditions. The main result of the present work is to identify initial conditions where 



2 Pratt's method to solve the hydrodynamical equations and to calculate the correlation functions differ 
in important aspects from our approach, but our final results for the correlation functions turn out to be 
rather similar to his. 
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the resulting stall is sufficiently long to be observable via pion interferometry |2T| , . As we 
show, the maximum time delay may occur at AGS as well as at RHIC energies, depending 
on whether matter is initially stopped in a fireball or expands as in Bjorken's longitudinally 
boost-invariant scenario. 

The remainder of this paper is organized as follows. In Section 2 we present the equation 
of state used in our investigations and explain our method to solve the hydrodynamical 
equations numerically for the above mentioned geometries. In Section 3 we compare the 
hydrodynamical expansion solutions for different initial energy densities, for different values 
of AT, and for different values of degrees of freedom in the QGP and hadronic phase. Section 
4 is devoted to the discussion of the lifetime of the system and how it can be measured via the 
widths of the side- and out-correlation function. Section 5 concludes with a summary and 
discussion of our results. Two Appendices contain the derivation of the formulae employed to 
calculate the two-particle correlation functions in the hydrodynamical framework. Natural 
units h = c = k B = 1 are used throughout this paper. 



2 Equation of State and Hydrodynamics 



Available lattice data for the entropy density in full QCD can be approximated by the simple 
parametrization [TJL I23J 



-CO 




tanh 



T-T r . 



AT 



where s c = const, x | (dq + dn) T c 3 is the entropy density at T c . Pressure p and energy 
density e follow then from thermodynamical relationships. For AT = 0, the equation of 
state (0) reduces to the MIT bag equation of state [^4| with bag constant B = | (dq/d H — 
1) T c s c / (dq/dn + 1)- If one measures energies in units of T c and energy densities in units 
of the enthalpy density e c + p c = T c s c , the equation of state (|l|) depends only on the ratio 
dci/d-H, and not on dq and dn separately. For AT = 0, this ratio determines the latent heat 
(density) e Q -e H = AB. (Here, e Q = |(4 d Q j d H -l)T c s c / (d Q / d H + 1) is the energy density 
at the phase boundary between mixed phase and QGP, e# = | T c s c / (dq/dn + 1) is that at 
the boundary between mixed and hadronic phase.) 

For the case dn — 3 (corresponding to an ultrarelativistic gas of pions) and dQ = 37 
(corresponding to two massless flavours of quarks and antiquarks, and eight massless gluons), 
the latent heat, 4B = 1.7 T c s c ~ 1.272 GeV fm -3 , is large. On the other hand, including 
a resonance gas in the hadronic phase and/or reducing the effective number of degrees of 
freedom on the QGP side [p5[| , dq/du = 3 may be taken as a (perhaps more realistic) lower 



limit, with a smaller latent heat €q — 6h = T c s c . Assuming that the high-temperature phase 



consists of gluons only (such as expected for the "hot-glue scenario" p6| ) this would then 
correspond to about 400 MeV fm~ 3 in physical units. To cover the range of uncertainty in 
the QCD equation of state we consider dq/dn = 37/3 and dq/dn = 3, and AT = and 
0.1 T c in our investigation of hydrodynamical expansion solutions. 

In Fig. 1 we show (a) the entropy density and (b) the energy density as functions of 
temperature, and (c) the pressure and (d) the velocity of sound squared c| = dp/de as 
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functions of energy density for AT = 0, 0.1 T c , and an ideal gas with dn degrees of freedom 
for dq/du = 37/3. The corresponding plot for dg/dn = 3 looks rather similar qualitatively 
and is not shown. Quantitative differences are: (i) the phase boundaries are different, 
€q shifts from 1.8125 T c s c to 1.375 T c s c , e# from 0.1125 T c s c to 0.375 T c s c , and p c from 
0.0375 T c s c to 0. 125 T c s c . Also, as mentioned above, (ii) the latent heat eg — tu is reduced. 

Figs. 1 (a,b) present the thermodynamic functions in a form to facilitate comparison 
with lattice data. Present lattice data for full QCD can be approximated with a choice 
of AT in the range < AT < 0.1 T c . In the hydrodynamical context, however, Figs. 1 
(c,d) are more relevant. As can be seen in (c), for AT = the pressure stays constant, 
p c = \T c s c / {dq/du + 1)) hi the mixed phase e# < e < €q. Hydrodynamical expansion 
is, however, driven by pressure gradients. It is therefore the velocity of sound, Fig. 1 (d), 
that is the most relevant measure of the system's tendency to expand. It represents the 
capability to perform mechanical work (which is proportional to pressure gradients dp) for 
a given gradient in energy density de. For AT = 0, the velocity of sound vanishes in the 
mixed phase, i.e., mixed phase matter does not expand at all on its own account, even if 
there are strong gradients in the energy density. This has the consequence that it does not 
perform mechanical work and therefore cools less rapidly. For finite AT, pressure gradients 
are finite, but still smaller than for an ideal gas equation of state, and therefore the system's 
tendency to expand is also reduced, cf. Fig. 1 (d). 

Hydrodynamics is defined by local energy-momentum conservation, 

dfjT^ 1 ' = . (2) 

Under the assumption of local thermodynamical equilibrium (the so-called "ideal fluid" 
approximation) the energy-momentum tensor T^ u assumes the particularly simple form |27) 

= (e + p) uV -pg^ , (3) 

where = 7(l,v) is the 4-velocity of the fluid (v is the 3-velocity, 7 = (1 — v 2 ) -1 / 2 , 
u^u^ 1 = 1), and g^ u = diag(+, — , — , — ) is the metric tensor. The system of equations (0) 
is closed by choosing an equation of state in the form p = p(e), i.e., as depicted in Fig. 1 
(c). In the ideal fluid approximation, the (equilibrium) equation of state is the only input to 
the hydrodynamical equations of motion @ that relates to properties of the matter under 
consideration and is thus able to influence the dynamical evolution of the system. The 
final results are uniquely determined once a particular initial condition and a decoupling 
("freeze-out") hypersurface are specified. 

The symmetry of the fireball and Bjorken cylinder geometry affects that the system of four 
equations (0) reduces to two independent equations. With the definition E = T 00 , M = T 0r , 
where the index r indicates the radial component of the energy-momentum tensor, the 
respective equations read 

d t E + d r [{E + p)v] = -F{E,p,v,r,t) , (4) 
d t M + d r (Mv+p) = -G(M,v } r } t) . (5) 

Here, v is the radial component of the velocity. For the fireball geometry, F and G do not 
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depend on the time t explicitly, 



F ib (E,p,v,r) = ^-{E + p) , G {h (M,v,r) = ^M . (6) 

For the Bjorken cylinder geometry, the above equations describe the system's transverse 
evolution at z = (cf. p8|), and due to the assumption of longitudinal boost invariance [20 



the hydrodynamical solution for arbitrary z can be easily obtained by a Lorentz boost with 
(space-time) rapidity rj = Artanh [z/t]. The functions F and G read in this case 

F m {E,p,v,r,t) = (H + iW + p) , G Bi {M,v,r,t)= ^ + iW . (7) 



In order to solve (^JS]), we employ Sod's operator splitting method p9"] , i.e., for each time 
step we first generate solutions of the hydrodynamical equations for F = G = 0. In this 
form, the equations are purely one-dimensional and can therefore be solved with e.g. the 
relativistic HLLE scheme presented in [|13|, The performance of this algorithm for solv- 
ing one-dimensional hydrodynamical problems and with equations of state featuring phase 



transitions was shown to be excellent [13], [3T| . It has also been employed in [14| to solve the 
Landau expansion problem. 

In a second step, Sod's method prescribes to correct E and M for the respective geometry 
by solving the ordinary differential equations 

^ = —F(E, p, v,r,t), ^ = —G(M, v,r,t) . (8) 

More specifically, in a slight variation of Sod's method we solve the finite difference equations 

E = E — F(E, p, v, r, t) dt , M = M — G(M, v, r, t) dt , (9) 

where quantities with a tilde are the solutions of the hydrodynamical equations with F = 
G = generated previously with the relativistic HLLE scheme. This two-step predictor- 
corrector scheme is repeated for each time step to obtain the complete time evolution of the 
system. 

Fig. 2 shows temperature and laboratory energy density profiles calculated with Sod's 
method for the expansion of (a,b) a sphere and (c,d) a Bjorken cylinder (for initial time 
to = r = 0.1 R) with an ideal gas equation of state p = e/3 in comparison to profiles 
generated with the semi-analytic method of characteristics E3] . That method is a benchmark 



test for numerical algorithms as long as the hydrodynamical solution is continuous |L3[ . The 
grid spacing for the HLLE scheme is taken as Ax = 0.01 R, the time step width for the HLLE 
scheme and the corrector step (|9|) is At = 0.99 Ax. As one observes, agreement is excellent 
(even on a logarithmic scale) which gives confidence that Sod's method as described above 
works for the more complicated equation of state (0) as well (which leads to discontinuous 
hydrodynamical expansion solutions for AT < AT* ~ 0.07676 T c (TT 
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3 Expansion solutions 



In this section we present the hydrodynamical expansion solutions for the fireball and Bjorken 
cylinder geometry for different initial energy densities €q. We assume eo to be homogeneous 
throughout the system. We compare solutions for the equation of state ([[]) with AT = and 
0.1 T c with solutions for an ideal gas of equation of state with du degrees of freedom. We 
explicitly show results for dqjdu = 37/3 and, where necessary, comment on the difference 
to the case dqjdu = 3. As in [jl4| we will use the notion "lifetime" for the intercept of a 
particular isotherm with the t— axis. 

Note that in our comparison we fix the initial energy density rather than the initial 
temperature or entropy density, as sometimes assumed |32| , because the latter are derived 
thermodynamic quantities while the former is determined by non-equilibrium energy loss 
mechanisms. For instance, at RHIC energies it is expected |26j that semi-hard perturbative 
QCD processes determine the initial energy density in the range eo ~ 10 — 20 GeV/fm 3 . 



3.1 Fireball geometry 

This case is rather similar to the one-dimensional Landau expansion studied in [Q. Dif- 
ferences are solely due to the spherical geometry (i.e., due to the extra terms (|6]) in the 
equations of motion). 

In Fig. 3 we show temperature profiles and isotherms in the t — r plane for 6q — e# = 
0.1125 T c s c . Of course, for this initial energy density, the case AT = (a,b) is identical to 
the ideal gas case (e,f). Note the delayed expansion in the case AT = 0.1 T c , Figs. 3 (c,d), 
due to the reduction in the velocity of sound, cf. Fig. 1 (d). This is quite similar to the 
one-dimensional expansion (cf. Figs. 3 (c,i) of Ref. [Q), although the overall scale of the 
lifetimes is now considerably reduced since the system has two more spatial dimensions into 
which it can expand. 

Fig. 4 shows the situation for an initial energy density eo = 1.875 T c s c which is close to 
tQ = 1.8125 T c s c . As in the one-dimensional expansion (cf. Fig. 5 of [|14[]), the lifetime of the 
system is very long in the case AT = 0, due to the small velocity of the deflagration front 
converting mixed phase matter into hadrons. As compared to the lifetime of the mixed phase 
in the one-dimensional expansion, tiif e ~ 40 i? (cf. Fig. 5 (g) of |14|), that lifetime is now, 
however, only about 18 R, due to the more rapid expansion in three spatial dimensions. As 
in the one-dimensional case, cooling is faster for AT = 0.1 T c (cf. Fig. 5 (i) of [0]). While 



in that case, however, the lifetimes were only slightly shorter than for AT = 0, now they 
are reduced by about a factor of 3. The fastest expansion is that for the ideal gas, in accord 
with the discussion of Fig. 1 (d). 

In Fig. 5 the initial energy density is well above the mixed phase region, eo = 18.75 T c s c ~ 
10 €q. As for the analogous one-dimensional expansion solution (Fig. 7 of [|14]]) the internal 



energy density and pressure are so large that, in the case of a transition to the QGP, Figs. 
5 (a-d), the system explodes rather than burns slowly as in the preceding case, Fig. 4. 
Nevertheless, the expansion is still somewhat delayed as compared to the expansion of an 
ideal gas, Figs. 5 (e,f). 
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3.2 Bjorken cylinder geometry 

For the discussion of the expansion in the Bjorken cylinder geometry it is instructive to first 
focus on purely longitudinal expansion. Longitudinal boost invariance with regard to boosts 
with n = Artanh [z/t] [^0] implies that the longitudinal velocity is given by v z = z/t. In 
turn, the longitudinal fluid rapidity is identical with the boost rapidity n. It is easy to show 
that in this case the hydrodynamical equations (§) reduce to 

T e|„ = -— , d v p\ T = 0, (10) 



where r = \Jt 2 — z 2 is the proper time associated with fluid elements moving with v z = z/t. 
The second equation implies that pressure gradients vanish along space-time hyperbolas 
t = const, and, for baryon-free matter, that the temperature is constant on these hyperbolas. 
The first equation is an ordinary differential equation on the space-time hyperbolas and 
describes the cooling of the system on account of the longitudinal motion. For baryon-free 
matter, the thermodynamical identities de = T ds and e + p = Ts imply the very simple 
cooling law 

ar «|, = ~, (ii) 

which has the solution s(r) = Sq [r/r ] _1 , independent from the underlying equation of state. 

For an ultrarelativistic ideal gas equation of state p — e/3, e = const, x T 4 , the solution 
of the longitudinal boost-invariant expansion problem is 

c(r) = e [t/to]- 4/3 , T(r) = T [t/t,]- 1 '* . (12) 

For the equation of state ([[]) with AT = 0, one can easily solve fllCf ) analytically. If eo > cq 
we have 

r (e - 5) [T/7-0]- 4 / 3 + 5 , r <r<r Q , 
e(» = S (eg + Pc) V/tq}- 1 - p c , r Q < r < t h , (13) 
[ e H [t/th]-^ , t>t h , 

f Tofr/ro]- 1 ^ , r < r < r Q , 
T(r) = T c , r Q < r < r H , (14) 

I T c [t/t h ]-V 3 , r>r^ . 

Here, rg = r [(eg — B)/ (eo — £>)]~ 3 / 4 is the (proper) time the system enters the mixed phase, 
and t h = TqdQ/dH the time corresponding to entry of the hadronic phase. Note that the 
time spent in the mixed phase is linear proportional to the ratio dqjdu- Also, the system 
does not cool at all in this phase, T = T c = const., due to the fact that no mechanical 
work is performed (cf. the above discussion of Fig. 1 (d)). As another consequence, the 
energy density does not decrease with r~ 4//3 as in the QGP and hadronic phase, but only 
proportional to r _1 , precisely as in a (one-dimensional) free-streaming expansion (where no 
mechanical work is performed as well). 
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For tQ > €q > tu , the solution reads 



c(r) = { (£o r + f e) 1 [T /? 1 " 1 " Pc ' T ° " T " TH ' (15) 
V ; \e H [t/t h ]- a / 3 , r > t h , v ; 




^(t) = , "-J ' (16) 



Here, r# = r (e +p c )/(^H +Pc)- For e < e#, the solution is identical to (|12"D . 

For finite AT there is, contrary to expectation, also a simple, semi- analytic solution 
which does not require to solve the equations of motion (|T0|) explicitly. Due to the fact 
that the entropy density behaves as s(r) = s [r/r ] _1 irrespective of the equation of state, 
and since the equation of state ([]]) establishes (for finite AT) a one-to-one correspondence 
between temperature and entropy density, there is also a one-to-one correspondence between 
temperature and r. Once s(r) and T(r) are known, one can easily calculate e(r) and p(r) 
from fundamental thermodynamical relationships. 

For illustrative purposes, we show in Fig. 6 the time evolution of (a) energy density, 
(b) entropy density, (c) pressure, and (d) temperature assuming an initial energy density 
eo = 10T c s c and oIq/cIh = 37/3. Solid lines are for AT = 0, dotted lines for AT = 0.1 T c , 
and dashed lines for an ideal gas with dn degrees of freedom. The change in the cooling 
law for the energy density when entering the mixed phase is clearly visible in Fig. 6 (a). 
Also, Fig. 6 (b) illustrates that the cooling law for the entropy density is independent from 
the equation of state, although the magnitude of the entropy density differs for the different 
equations of state because of our assumption of fixed initial energy density. 

We mention that for AT = 0.1 T c , the parts of the system with temperature T > 0.8 T c 
cool faster than for the ideal gas, cf. Fig. 6 (d). Only parts cooler than ~ 0.7 T c survive 



longer than in the ideal gas case. This is in agreement with the results of |1J] for the Landau 



expansion model, where it was concluded that it is therefore less likely that electromagnetic 



radiation is a viable signature for the transition, in contrast to the conclusions of On 
the other hand, it was speculated in |nj (referring to the conclusions of [fZT], that, if the 
system freezes out at temperatures < 0.7 T c , the delayed expansion might be observable in 
side- and outwards radii of two-particle correlation functions. This indeed happens to be 
the will be shown in the next section. 

For given initial energy density eo and proper time r , one can ask the question, how 
much longer it takes a system described by the equation of state (|l|) to reach a certain 



freeze-out temperature Tf < T c in comparison to the expansion of an ideal gas (12). Since 



the ideal gas is supposed to freeze out at the same temperature T/ and to have da massless 
degrees of freedom, particle and entropy density at freeze-out are the same in both cases. 
For AT = one can easily answer this question analytically: for eo 3> B, the freeze-out time 
Tf = r(Tf) is prolonged by a factor [c^q/g^] 1 / 4 . For dQ/d H = 37/3, the delay in cooling is 
about 87% (as one can also see in Fig. 6 (d)), for dq/du = 3 it is only about 32%. Note that 
this prolongation of the lifetime of the system is in agreement with the one for the Landau 
expansion model at very high initial energy densities (cf. JT3] and Fig. 10 below). Only for 
initial energy densities around €q is the prolongation of the lifetime larger in the Landau 
expansion (see Fig. 10) due to the fact that there exists a slow deflagration solution which 
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delays the expansion even further. Such a solution does not occur in the purely longitudinal 
Bjorken expansion. 

This changes when we consider transverse motion as well. In the following we fix tq = 
0.1 R. This choice is motivated by the fact that for Au+Au collisions at RHIC, the transverse 
radius of the hot zone is of the order 5 fm, while the time scale of local equilibration is roughly 
given by the energy loss of a parton in strongly interacting matter, r dE / dx ~ 0.5 fm |TT|. We 



have, however, also considered the case To = 0.5 R and r = 1/3 T (which is motivated by 
an uncertainty principle argument) and will comment on differences to the choice r = 0.1 R 
where necessary. 

Let us first consider the case eo = eu = 0.1125 T c s c , cf. Fig. 7. Due to the strong initial 
longitudinal motion, the system cools rather quickly below temperatures of 0.5 T c , even before 
the transverse rarefaction wave reaches the center of the cylinder. The lifetimes are therefore 
solely determined by the longitudinal scaling expansion. This causes the horizontal parts 
of the isotherms. For the T = 0.5 T c isotherm and AT = one expects from eq. (0) the 
lifetime tnf c = 0.5~ 3 x 0.1R = 0.8 R, in good agreement with Figs. 7 (a,b,e,f). (Due to the 
numerical dissipation adherent to any finite difference scheme that solves the hydrodynamical 
equations, the system cools slightly faster in the numerical calculation, cf. Figs. 7 (b,f).) The 
effect of longitudinal cooling is reduced for larger values of To, cf. (0). For instance, in the 
case To = 0.5 R the center still has T ~ 0.6 T c when the transverse wave reaches it. The 
lifetimes for T > 0.6 T c are therefore (as expected) a factor of 5 longer. 

In Fig. 8 we show the situation for eo = 1.875 T c s c ~ eg. As in the preceding case, it 
takes some time before the transverse expansion wave reaches the center. Thus, the lifetimes 
of the highest temperatures in the system are solely determined by the longitudinal scaling 
solution presented above (horizontal parts of the respective isotherms). 

It is interesting to note that the lifetimes for AT = and this particular eo are not 
that exceedingly long as they were in the case of one-dimensional and spherical expansion. 
These long lifetimes were due to the slow velocity of the deflagration front. A deflagration 
solution, however, exists only for energy densities in the mixed phase, where matter is 
thermodynamically anomalous ||14|| . Here, the strong longitudinal motion cools the system 
quickly below the respective energy densities, such that the transverse expansion proceeds as 
a (comparatively fast) simple rarefaction wave instead of the deflagration. Thus, the system 
cools even more quickly and the lifetimes are considerably reduced. 

This effect is, however, compensated when higher initial energy densities are considered. 
In Fig. 9 we show the solutions for eo = 18.75 T c s c ~ 10 cq. This case is close to the initial 
conditions expected from mini-jet production at RHIC energies. Now the initial energy 
density is sufficiently high that for AT = the plasma enters the mixed phase at relatively 
late time (tq ~ 3 fm for r = 0.5 fm, see also Fig. 6 (d)). Therefore, the time spent in 
the mixed phase is long enough in this case to allow a transverse deflagration front to form, 
cf. Figs. 9 (a,b). In comparison to the ideal g (e,f), where no such solution exists, 

the expansion is therefore considerably prolonged. Such a time delay is observed also for 
AT = 0.1 T c , Figs. 9 (c,d), and is here due to the reduction of the velocity of sound in the 
transition region. 

We note that for tq = 0.5 R, the initial energy density where one observes this effect is 
smaller, eo ~ 2 — 3 eg. On the other hand, for tq = 1/3 Tq the initial time decreases with 
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increasing initial energy density, enhancing the effect of longitudinal cooling, such that eo has 
to be very high (> 40 eg) for the onset of this effect. Then, however, a steady increase of the 
initial energy density affects a steady decrease of tq such that one observes a prolongation of 
the lifetimes over a wider range of energy densities (see also Fig. 13 below). It is furthermore 
interesting to note that this phenomenon occurs at about the same initial energy densities 
(in units of T c s c ) in the case dg/du = 3 (cf. Fig. 12). 

Finally, for very high energy densities, one recovers the case of the violent explosion, 
and the lifetimes are again reduced. This happens at lower initial energy densities for the 
case tq = 0.5 R, since the longitudinal motion is less effective in cooling the system. For 
To = 1/3 To we have not found a decrease in the lifetime up to the highest initial energy 
density studied by us, eo = 300 T c s c . The reason is, as mentioned above, that the increasingly 
more efficient cooling due to longitudinal motion compensates the increasing tendency of the 
system to explode transversally. 



4 Lifetimes and two— particle correlations 

In this section we first discuss the lifetimes of differently hot parts of the system (defined 
as the intercept of a particular isotherm with the t-axis) as a function of the initial energy 
density. Assuming that the system freezes out at these temperatures, we then show how the 
lifetimes can be inferred from two-particle correlation functions. 



4.1 Lifetimes 

For convenience, we first show in Fig. 10 the lifetimes as a function of eo for the one- 
dimensional (Landau) expansion studied in |T4J. Figs. 10 (a) and (b) are for the case dg/dn 



37/3, (c) and (d) are for dg/du = 3. The thick lines in (a,c) show the lifetimes of matter 
with T = 0.7 T c (solid line), T = 0.9 T c (dotted), and T = T c (dashed) for AT = 0, the 
corresponding lines in (b,d) are for finite AT = 0.1 T c . For comparison, the thin lines are the 
corresponding lifetimes for the expansion of an ideal gas with dn degrees of freedom (they 
are identical in (a) and (b) as well as in (c) and (d)). 

One clearly observes the maximum in the lifetime of mixed phase matter (T = T c ) at 



eo = eg in Figs. 10 (a,c). As discussed in |TJ]], this maximum vanishes for finite AT, Figs. 



10 (b,d). Moreover, such hot matter does not survive as long as in the ideal gas case. The 
lifetime of cooler matter with T = 0.7 T c , however, stays longer than in the expansion of an 
ideal gas, independent of the value of AT. Note the change in the scale of tm e by a factor of 
3 between (a,b) and (c,d). This is due to the reduction of the latent heat for a smaller ratio 
dg/dn, which in turn accelerates the expansion. 

Fig. 11 shows the corresponding diagram for the spherical expansion. First of all, one 
notices that, at high eo, there is almost no difference between the lifetime of a system with a 
transition to the QGP and an ideal gas. Inspecting Fig. 5, one observes, however, that the 
bulk of matter at finite r does indeed live longer in the first case. Thus, for the spherical 
expansion, our definition of 'lifetime' should rather be replaced by an average over the 
particular isotherm. It is, however, not necessary to do so at this point, since the correlation 
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functions considered subsequently will take this into account in a natural way. 

Second, one observes that the lifetimes do not grow as strongly for high eo as in the one- 
dimensional expansion. This is due to the fact that the system disperses its initial internal 
energy much more efficiently into kinetic energy in three dimensions than it does in one 
dimension. This also leads to the reduction of the overall scale in the lifetimes as compared 
to Fig. 10. 

Third, it is noticeable that the increase of the lifetime for AT = at eo = eg (where the 
lifetime is maximum) as compared to the ideal gas case is about a factor of 2 bigger than in 
the one-dimensional expansion. The prolongation of the lifetime is thus most pronounced 
in spherical geometry. Unfortunately, the reduction of the lifetime in the case of a smooth 
transition is also rather strong, cf. Figs. 11 (b,d). Nevertheless, for freeze-out at T = 0.7 T c 
the lifetime can still be longer by a factor of 2 in the case of a transition than in the expansion 
of an ideal gas. Moreover, in contrast to the one-dimensional case, the spherical geometry 
leads to a (broad) maximum in the lifetime around €q. 

In Fig. 12 we present the lifetimes for the expansion of a Bjorken cylinder with initial 
time r = 0.1 R. One still observes the distinguished maximum in the lifetime associated 
with the transition to the QGP as seen in the preceding Fig. 11. In this case, however, the 
overall scale is even smaller (due to the fact that the system is from the very beginning quite 
effectively diluted by the longitudinal velocity field). 

Moreover, one observes a shift in the maximum of the lifetime. This was discussed in 
the preceding section and is due to the fact that the longitudinal motion compensates the 
tendency to explode transversally (on account of a high initial energy density). For AT = 
the maximum is now around eo ~ 20 — 50T c s c for both values of oIq/oIh- For AT = 0.1 T c , 
one still observes a maximum in the lifetimes, but as in Fig. 11 it is also broader and less 
pronounced. 

The corresponding plot for To = 0.5 R looks rather similar, and will therefore not be 
shown here. The position of the maximum is, however, shifted to smaller eo ~ 3 — 7T c s c , 
for reasons discussed above. The situation for a dynamical tq = 1/3 To is shown in Fig. 13. 
As one expects from the above discussion, there is a steady increase of the lifetime up to the 
highest considered e . 



4.2 Two— particle correlation functions 



To calculate two-particle correlation functions we use the method developped by Pratt [^T 



Sinyukov [[33] and others, and applied to hydrodynamics by the Marburg group |3"3 |. This 



method is essentially a straightforward generalization of the Cooper-Frye formula |35| for 
single inclusive particle spectra to the case of correlation functions. 

To calculate single inclusive particle spectra or two-particle correlation functions in the 
hydrodynamical framework, one assumes that a fluid element decouples from the fluid evolu- 
tion ( "freezes out" ) as soon as the particle density drops below a certain critical value where 
the collision rate becomes too small to maintain (local) thermodynamical equilibrium. In our 
case, the particle density depends on the temperature T only. Therefore, freeze-out happens 
across isotherms in the space-time diagram. To incorporate the freeze-out consistently into 
the solution of the hydrodynamical equations is a non-trivial problem. Hydrodynamics is 
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obviously not applicable to determine the motion of particles that are already frozen out. 
Therefore, the hydrodynamical solution has to be restricted to a limited region of space- 
time. The main difficulty is that its boundary, i.e., the freeze-out hypersurface, has to be 
determined dynamically, i.e., simultaneously with the solution of the hydrodynamical equa- 
tions. A possible treatment of this problem was recently proposed by Bugaev [BJJ, but has 
so far not been applied in practical calculations. 

Commonly one circumvents this problem assuming the validity of the hydrodynamical 
description in the whole forward light cone, then solves the hydrodynamical equations, and 
determines the isotherms. Finally, one employs the Cooper-Frye formula ]35] to calculate 
particle spectra along the isotherm corresponding to the particular freeze-out temperature 
(for a detailed discussion of this approach see, for instance, [0). We note that the Cooper- 
Frye formalism holds rigorously only on space-like hypersurfaces. For time-like hypersurfaces 
it may yield negative numbers of particles emitted from the hypersurface (corresponding to 
particles that do not freeze out but reenter the fluid). To cure this obviously unphysical 
result, modifications have been proposed in [BJ], j3"5| . The approach of Ref. |38| is, how- 
ever, problematic since it does not describe correctly radiation from a static source. When 
modifying the particle spectra as proposed in [BJJ we found, however, for the cases under 
consideration virtually no deviation to the results obtained with the conventional approach 
of Cooper and Frye. Therefore, we will use this well-established method throughout the 
following. 

The two-particle correlation function measures the coincidence probability P(pi,p2) 
of two (identical) particles with momenta pi, p 2 relative to the probability of detecting 
uncorrelated particles from different events, 



C2(Pl,P 2 ) 



^(Pl,P2) 



(17) 

(Pi +P2V2 and the relative 



Pi ~~ P2- Under the assumption that the particle source is chaotic and 



In the following, the average 4-momentum is denoted as 
4-momentum as 

sufficiently large, and that the emitted particles are bosons (with degeneracy factor d) the 
two-particle correlation function can be written as [Bl] 



C2(Pl,P 2 ) = 1 + 



(27T)- 



J E dS-K exp[i£.g] 



dN 



Pi 



IF 



rr dN 



where [Bo" 



E 



dN 



d 



d 3 p (2vr) 3 ^ J \ T ) (19) 

is the single inclusive momentum distribution. f(x) = (e x — l)" 1 is the Bose-Einstein 
distribution function, and the fluid 4-velocity. The integrals run over the freeze-out 
hypersurface. In general, that hypersurface is represented by a 3-parametric (4-vector) 
function £ M (£, rj, <fi), and the normal vector on the hypersurface is determined by 



d£, 



«9£ Q <9£t 
d( drj dtp 



d( dr]< 



(20) 
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where e lxa p 1 (= —1 for (fia^j) an even permutation of (012 3)) is the completely antisym- 
metric 4-tensor. While the (freeze-out) temperature T is (by definition) assumed to be 
constant along the freeze-out hypersurface, the fluid velocity varies, w M = m^(S). 

For symmetric systems, the number of independent variables which C 2 depends on can 
be reduced. For the spherically symmetric fireball geometry, the orientation of the average 
momentum can be chosen arbitrarily, such that C 2 depends only on the modulus of K. 
Our choice of coordinate system will be such that K = (0,0,7^). Furthermore, q can be 
decomposed into a so-called "out" component q out = (0, 0, q Q ut) parallel to K and a "side" 
component q S id e orthogonal to K and q ou t- Rotational symmetry around the direction of K 
allows us to choose q S id e = (<? S ide, 0, 0), such that the correlation function depends only on 
three independent variables, C 2 (K, g out , g S ide)- 

For the Bjorken cylinder geometry, we restrict our consideration to particles emitted at 
midrapidity, K z = q z = 0. Rotational symmetry around the 2-axis allows us to choose 
the average transverse momentum as Kj_ = (K, 0,0), and consequently, q out = (g out ,0,0), 
q S ide = (0, q S ide, 0). Again, C 2 (K, q OVLt , q S id e ) is a function of three independent variables only. 
The explicit evaluation of for the fireball and Bjorken cylinder geometry is referred to 
Appendices A and B, respectively. 

As is well known |2]J, the width of the correlation function in out-direction is inversely 
proportional to the duration of particle emission, i.e., to the lifetime of the source. Anal- 
ogously, the inverse width of the correlation function in side-direction is a measure for the 
(transverse) size of the source. To be more precise, for fixed average (transverse) momen- 
tum K we define side- and out-correlation functions as C 2 , side (? side) = C 2 (K, 0, g S ide) and 
C2,out(<?out) = C 2 (K, q ont , 0), respectively. We furthermore define a corresponding inverse 
width as -R S id = V^sido where g* ide is determined by C 2 ,sidc(<?side) = 1-5, and analogously 
for ^0^. 

We emphasize that we do not attempt a standard Gaussian fit of the two-particle corre- 
lation function, which would relate the inverse widths to the usual radii parameters. First 
of all, we will be interested only in the generic shape of the correlation functions which is 
satisfactorily characterized by their inverse widths -R S id e , -Rout- Second, the functional form 
of the correlation functions is not a Gaussian (cf. the explicit formulae in Appendix A and 
B). The fit procedure would therefore only introduce unnecessary errors. Third, the radii 
parameters as well as the inverse widths are only proportional to the actual (average) size 
and lifetime of the system. It is well known that hydro dynamical flow affects the radii pa- 



rameters/widths, and, moreover, that this effect is sensitive to the choice of K [21 , |34|| . It 



is therefore tedious (if not impossible) to relate -R S id e and _R out to the real source size and 
lifetime. This holds for our model calculations as well as for the experiment. 

We can expect, however, that such effects are either irrelevant or largely cancel out if 
we consider the ratio _R ut/-Rside- As can be inferred from the space-time diagrams in Figs. 
3-5 and 7-9, while the size of the system is approximately constant, it is the lifetime that 
varies appreciably, depending on whether the system undergoes a phase transition or not. 
Thus, -R ou t/-Rside can be expected to be a good measure for the lifetime of the system. In 
the following, we choose the particles to be pions with mass m = 138 MeV, and as average 
momentum K = 300 MeV. On one hand, this value is well within the typical experimental 
acceptance. On the other hand, it is large enough to exhibit effects of the prolonged lifetime 
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in the case of a transition to the QGP pi| . To fix the g-scale in MeV we take R = 5 fm as 
a typical initial radius of the system. 

In Fig. 14 we show as an example side- and out-correlation functions corresponding to 
the hydrodynamical evolution of Fig. 4. In light of the uncertainties in the actual freeze-out 
temperature, we calculate the correlation functions along isotherms for T = T c , 0.9 T c , as 
well as 0.7 T c . First of all, one immediately recognizes that in all cases the widths of the 
correlation functions correspond closely to the space-time structure of the corresponding 
isotherms. For instance, the exceedingly long lifetime of the system in Fig. 4 (b) is reflected 
in the comparatively small width of the out-correlation function Fig. 14 (b). On the other 
hand, the transverse size is about the same in all cases, cf. Figs. 4 (b,d,f), which reflects 
in nearly identical side-correlation functions in Figs. 14 (a,c,e). Furthermore, with the 
exception of the case AT = 0.1 T c , the correlation functions are almost completely insensitive 
to the freeze-out temperature chosen. This is intuitively clear since the corresponding space- 
time isotherms differ only marginally for AT = and the ideal gas calculation, cf. Figs. 
4 (b,f), while there are larger differences in Fig. 4 (d). Note that the results shown in 
Figs. 14 (a,b,e,f) correspond nicely to those in Fig. 4 of [^T|, in spite of the fact that the 
latter calculation employs inhomogeneous initial conditions, a different way to solve the 
hydrodynamical equations, and a different treatment of the freeze-out^} In Fig. 15 we show 
the correlation functions corresponding to the hydrodynamical evolution of Fig. 9. Again, 
they adequately characterize the space-time geometry of the source. 

Fig. 16 shows the experimentally measurable ratio R out / Rside as a function of eo for the 
spherical fireball geometry. Comparing the results with Fig. 11, one observes that this ratio 
reflects closely the behaviour of the lifetime of the system, independent of details in the 
equation of state such as the width of the transition region AT or the latent heat of the 
transition (which is proportional to dg/dn)- Also, for the case of a first order transition, 
AT = 0, Figs. 16 (a,c), the enhancement in R out /R S ide over the ideal gas case is a factor of 
3 to 7 (for dq/dn = 3 to 37/3) at e ~ €q. In the case of a smooth transition, AT = 0.1 T c , 
Figs. 16 (b,d), this is considerably reduced (as expected from Fig. 11), but if the system 
freezes out at temperatures Tf < 0.7 T c , there is still a factor of 2 enhancement over the 
ideal gas case. 

In Figs. 17, 18 we present the corresponding results for the Bjorken cylinder expansion 
with t = 0.1 R and r = 1/3 T , respectively. In all cases we find that the experimentally 
measurable ratio of correlation widths mirrors closely the dependence of the lifetime on 
initial conditions in Figs. 12, 13. The most favourable case R ou t/R S idc ~ 3.5 (for a strong 
first order transition with a large latent heat, dqjdu = 37/3, cf. Fig. 17 (a)), may be reached 
with initial conditions expected at RHIC energies. At that point the enhancement of this 
ratio is about a factor of 2 above the ideal gas case and virtually independent of the freeze- 
out temperature. In the scenario with tq = 1/3 To, Fig. 18 (a), the enhancement is somewhat 
smaller and is shifted toward higher initial energy densities because the initial pure Bjorken 
expansion phase starts earlier than in the fixed tq = 0.1 R case. For a smooth transition with 
AT = 0.1 T c , the maximum ratio is reduced to 2.5 and varies less rapidly with eo, but is still 

3 The freeze-out in |2lJ is in some sense performed dynamically, the freeze-out surface, however, is assumed 
to have no time-like parts. 
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about 40% larger than for the ideal gas expansion. In this case, however, a significant time 
delay can only be observed if the freeze-out occurs relatively late with Tf ~ 0.7 T c . As seen 
in Figs. 17, 18, for earlier freeze-out, smaller To, or smaller dq/dn the enhancement relative 
to the ideal gas case is significantly reduced and would be more difficult to observe. 

At energy densities estimated to be reached in CERN SPS Pb+Pb-collisions (e ~ 1 — 
2T c s c in our units), we expect from our results that -R ou t/-Rside ~ 1-5 — 2. However, present 
data from CERN SPS |39| indicate that the (fitted) out-radii are rather similar to the side- 
radii. This does not contradict our results, because, as shown by Schlei et al. HO] in the 
framework of a hydrodynamical calculation similar to ours, correlation functions constructed 
from thermal pions only give i? out /i? S id e ~ 2 (cf. especially [j34}Q, while the incorporation 
of long-lived resonance decays leads to a reduction of that ratio and good agreement with 
the measured radii. We note that kaon interferometry f4l| , EH, EH1 * s preferable, though 
experimentally more difficult, because only distortions of the interference pattern due to 
shorter lived K* resonances have to be taken into account. 



5 Conclusions 

In this paper we investigated the spherically symmetric expansion of a fireball as well as the 
cylindrically symmetric transverse expansion of a QGP with boost-invariant initial conditions 
along the beam axis. The expansion was treated in the framework of ideal relativistic 
hydrodynamics and extends our systematic study [13. 14. 16. 31, 37] of collective flow patterns 



with realistic equations of state. The symmetries of the considered geometries allowed us 
to calculate the flow patterns using a simple modification of a well-tested one-dimensional 
algorithm. 

The emphasis of the present investigation was on how a rapid cross-over to the QGP 
in the equation of state influences the collective expansion dynamics. In particular, since 
present lattice data only constrain the width of the transition region to the QGP to be in 
the range < AT < 0.1 T c , it is important to test how such uncertainties may influence 
dynamical observables. We also studied the effect of varying the ratio of degrees of freedom 
in the QGP and hadron phase (i.e. essentially the latent heat of the transition) on the 
system dynamics. The results were compared to the expansion of an ultrarelativistic ideal 
gas (without transition) using fixed energy density initial conditions. 

We focussed on the lifetime of the system as a function of initial energy density as an 
important collective observable that can discriminate between different equations of state. As 
expected from previous one-dimensional studies in the framework of the Landau expansion 



model ||14|| , we found that the lifetime of a spherical fireball is much longer in the case of a 
first order phase transition, AT = 0, as compared to the expansion of an ideal gas without 
transition. The prolongation of the lifetime in that case can be up to a factor of 4.5 to 9 (for 
dqjdu = 3 to 37/3, respectively), provided the initial energy density corresponds to that of 
mixed phase with a large fraction of QGP. This time delay effect was originally pointed out 
by Pratt [^TJ. In the case of a smooth transition, AT = 0.1 T c , however, this time delay is 
drastically reduced. The lifetimes (of matter with T = 0.7 T c ) nevertheless still remain on 
the order of a factor of 2 longer as compared to the ideal gas expansion. 
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For the Bjorken cylinder expansion, it is also necessary to specify the initial (proper) time 
To for the onset of hydro dynamic expansion. To explore uncertainties associated with this 
additional degree of freedom, we investigated the cases r = 0.5 R, 0.1 R, and a dynamical 
t = 1/3 To varying with the initial temperature. The results were similar as for the spherical 
expansion, up to two important exceptions: (a) the maximum lifetimes emerged at higher 
initial energy densities (the exact value of which depends on the choice of To) corresponding 
to QGP matter instead of mixed phase matter, and (b) the lifetimes were in general shorter. 

Both effects are explained by the very efficient cooling due to the initial longitudinal 
velocity profile associated with the boost invariance of the problem. This effect causes an 
overall reduction of the lifetimes. Moreover, in order for slow (van Hove deflagration 
fronts to dominate the cooling mechanism in transverse direction, one has to start at higher 
initial energy densities to compensate for the longitudinal cooling. Otherwise, the longitu- 
dinal cooling reduces the energy densities too fast and the associated deflagration solution 
13fl vanishes. Further cooling is then achieved (in the usual way) through a (fast) simple 



wave. 

Finally, we showed (cf. also pi, 22) that the prolongation of the lifetime can be ob- 



servable via the ratio R ou t/ R s ide of inverse widths of two-particle correlation functions in 
out- and side-direction. This ratio follows the behaviour of the lifetimes rather closely. The 
prolongation of the lifetime in the case of a transition to the QGP could therefore be in 
principle searched for using this observable. The enhancement of that ratio is, of course, 
strongest in the case that the transition is first order with a large latent heat. An interesting 
result is that, for the fireball geometry, the effect is maximum for energy densities achieved 
at the AGS, while for the Bjorken cylinder geometry, the maximum of R ou t /R side occurs at 
energy densities presumably reached at the RHIC collider. 

There are several effects which may reduce the strength of the time-delay signal observ- 
able via the i? out /-R S ide~ r atios that will require further investigation. First, for the fireball 
geometry it is important to extend our considerations to finite baryon number density. At 
finite chemical potential the width of the cross-over region may be significantly larger than 
at zero chemical potential. Second, the decay of long-lived resonances can simulate time 
delay f44 |. Interferometry with kaons instead of pions is therefore preferable |43|. Finally, 
while our investigations covered a wide range of uncertainties in the equation of state, our 
calculations have neglected effects of dissipation that tend in general to reduce the collective 
flow strengths predicted via ideal hydrodynamics. For instance, bulk viscosity appears in the 
hydrodynamical equations of motion in a similar way as the pressure, and could in principle 
counteract any reduction of the velocity of sound in the transition region. The main result 
of this paper is, nevertheless, that the generic time-delay signature of QGP formation is 
remarkable robust to present uncertainties in the QCD equation of state. 
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Appendix A 



In this Appendix we explicitly calculate the out- and side-correlation functions in the fire- 
ball geometry. We first specify the parameters of the freeze-out hypersurface. For the 
following discussion, we conveniently choose spherical coordinates r, 9, <p in space. Let (a) ( 
parametrize the surface in the t — r plane, < ( < 1, with the lower boundary corresponding 
to t = and the upper to r = 0. Let (b) rj parametrize the surface in the r — 9 plane and (c) 
in the r — ip plane. In the t — r plane, the surface is given by the corresponding isotherm 
shown in Figs. 3-5. Due to spherical symmetry, however, the surface in the r — 9 and r — ip 
planes is trivial. We can simply identify 7] = 9, <p = ip and vary them within the standard 
boundaries for the azimuthal and polar angle. As a consequence, the time tf and the radius 
rf at freeze-out depend only on (, not on rj or <p. The freeze-out surface is thus given by 
= (i/(e)j ff{0 e r ), where e r = (sinr? cos 0, sin 77 sin 0, cos?]). 

Applying (^) and choosing the sign to have dE M pointing outwards along the isotherm, 
we obtain 

dS M = ^-^,^e r ) rJ(C) sin?? dCd??d0 . (21) 

For the calculation of the correlation function we employ the Boltzmann approximation. 
For an average momentum K = 300 MeV chosen in our calculations the error introduced is 
negligible, in particular since we also consider the particles to be massive pions, m = 138 
MeV, and choose T c = 160 MeV. 

For the single inclusive pion spectrum (|T9|), one conveniently chooses p = (0,0,p), and 
obtains after inserting (|2l| ) and performing the 0- and ^-integrations 



_, diV d r 1 n> 9/A . p < T f sinha drv 



cosh a sinh a 



a 2 



(22) 
-1/2 



where a = pvy/T, E = [p 2 + m 2 ] 1 / 2 , and v is the (radial) fluid 3- velocity, 7 = [1 — 
Here (and in the following), the ^-integration has to be done numerically along the respective 
isotherms obtained from the solution of the hydrodynamical equations. For further use, we 
define 

For the side-correlation function, g out = 0, and consequently 

pt = (£,g sidc /2,0,/0 , (24) 

f 2 = (E,-q side /2,0,K) , (25) 

K» = (£7,0,0,*) , (26) 

<? = (0,g sidc ,0,0) , (27) 

where E = [K 2 + g 2 ide /4 + m 2 ] 1 / 2 . Since the single inclusive spectrum ( p2[) does not depend 
on the direction of p, the two single inclusive spectra in the denominator in (|T^) are equal. 
Furthermore, for the numerator one has to calculate the expressions 

Ji = / dS • Ke~ u - K/T cos [E • q] , (28) 
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T 2 = J^-K e - u - K/T sin [£•<?] . 



(29) 

For the calculation of X\ we insert (j2~I|) and perform the 0-integration with the help of 
eq. (3.715.18) of Ref. |45l| . Subsequently, the ^-integration can be done with an analytic 
continuation of either eq. (6.616.5) or (6.677.6) of and a suitable first derivative of these 
formulae. The final result reads 
-l 



2i 



47T 



dCr^(C) e 



-S7/T 



(30) 



sinh \J a 2 — b 2 drf 

~~ P d^ + v 7 ^ 



v 7 ^ 



6 2 



cosh y/ a 2 — b 2 sinh \/a 2 — b 2 



Va~ 2 



6 2 



6 2 



where a ee Kv^f/T, b = q S ide^f- It is easy to show that for the side-correlation function 
1 2 ee by symmetry. The final result thus reads 

C*2,side — 1 +2^l/^0 J 



with 1q from ( p3|) and Xi from ( ]30| ) . 

For the out-correlation function, g side 



Pi 



0, and consequently 

(#1,0,0,^ + ^/2) , 
(£ 2 ,0,0,K-g out /2) , 
(AT , 0,0, AT) , 
(E 1 - E 2 , 0, 0, q ont ) , 



(31) 



(32) 
(33) 
(34) 
(35) 



where E 1)2 = [{K±q out /2) 2 + m 2 } 1 ^ 2 , K° = (E 1 + E 2 )/2. The 0-integration in the calculation 
of X\, 1 2 is now trivial. For the ^-integration, we employ the angle addition theorem and 



eqs. (2.662.1,2) of @], obtaining 

Ji = 47r^ 1 dCrJ(C)e-^° 7/T |-K [cos /3 Jo(a, 6) + wa/3J 1 {a i b)] 



+AT 

where (3 = (Ei — E 2 )tf, and 

a cos 6 sinh a + b sin 6 cosh a 
Jo{a,b) ee 



o dJ (a,b) . R dJ x {a,b) 
cos p h sin p 



da 



da 



dtf 
dC 



~TF , (36) 



Ji(a,6) 



a sin 6 cosh a — b cos 6 sinh a 



a 2 + 6 2 ' n " a 2 + 6 2 

with a ee Kv~f/T and 6 ee g out r^. In complete analogy one obtains Z 2 as 

X 2 = 47r^ 1 dCrJ(C)e-^/ T |-K° [sin (3 J (a,b) - cos (3 J^b)} ^ 



(37) 



. dJ Q (a,b) R dJ x {a,b) 
sin p - cos p 



da 



da 



dC 



Finally, the out-correlation function reads 

C2,out = 1 + 



J (Pl) ^o(P2 



(3* 



(39) 



20 



Appendix B 



In this Appendix we explicitly calculate the out- and side-correlation functions in the 
Bjorken cylinder geometry. Convenient coordinates to work in are cylindrical coordinates 
r, if, z in space. As in the preceding Appendix, (a) £ parametrizes the freeze-out surface in 
the t — r plane (at z = 0), < ( < 1, t/(0) = r , r/(l) = 0. (b) r\ parametrizes the surface 
in the t — z plane. Due to boost invariance, that surface is simply the space-time hyperbola 



of constant r/ = Jt^ — zj. Note that at z = 0, Tf = t/(C) an d therefore Tf depends on (. 
A natural choice for the parameter 77 is the space-time rapidity Artanh [zf/tf] (as implied 
by our choice of symbols), such that tf((,rj) = Tf(Q cosh?], Zf((,r)) = Tf(Q sinhr?. Finally, 
(c) parametrizes the hypersurface in the r — if plane and due to cylindrical symmetry can 
be identified with the polar angle, = ip. We note that due to boost invariance along z, Tf 
cannot depend on 77. 

The hypersurface 4-vector thus reads 

SM = ( T /(C) cosh??, r/(C) cos0, 77(C) sin0, r/(C) sinh??) , (40) 

with T/(C) = tf((, 77 = 0), i.e., given by the isotherms in Figs. 7-9. The normal vector (with 
proper orientation) is readily calculated with 



dS^ = f-^r cosh??, ^ cos0, ^ sin0, ^ sinh 77 J r f (()T f (() d(dr]d(f) . (41) 

For the calculation of the single inclusive spectrum, we employ m± = [pj_ + m 2 ] 1 / 2 , i.e., 
p° = mj_ cosh 7/, = mj_ sinhy, with the longitudinal (particle) rapidity y = Artanh [p z /i?]. 
Furthermore, due to rotational symmetry around z, we may choose = (p±, 0, 0). Thus, 

dS ■ p = ^-m_L cosh [y - 77] ^ +p±cos<f) 77(C) r/(C) dCdr?d0 . (42) 

The fluid 4-velocity in the Bjorken cylinder expansion reads = 7(1, v e r , z/t), where 
v = vj_, e r = (cos 0, sin 0, 0). With the (longitudinal) space-time rapidity 77 = Artanh [z/t] 
and the transverse fluid rapidity r] r = Artanh [v cosh?]], one obtains 

= (coshr? cosh?/ r , sinh7/ve r , sinh 77 cosh/7 r ) . (43) 

Note that longitudinal boost invariance implies that r] r cannot depend on 77, and is therefore 
given by the solution of the hydrodynamical equations at z = 0, r] r = Artanh [v(z = 0)], as 
provided in Section 3. 

For the single inclusive momentum distribution (in Boltzmann approximation) we insert 
(|^) and fl43|) into (|i~9f), perform the //-integration with the help of eq. (3.547.4) of f|5j, and 



the (^-integration using the formula (3.937.2) |]4"5fl , resulting in 

'p± sinh rj r \ drf 



+p±Ko Jru^tj 7j /Pxsmh^A M (44) 
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Note that the final spectrum respects the symmetries of the problem, i.e., it is azimuthally 
symmetric and boost invariant along z. For further use, let us define 



2b (P j 



(2vr) 3 E dN 



d 



d 3 p 



(45) 



For the side-correlation function for particles at y = we may choose = (K, 0, 0) such 
that 



Pi 
u 
P2 



(E, K, g side /2, 0) , 
(E, K, -g sidc /2, 0) 
(E,K,0,0) , 
(0, 0, g side , 0) , 



(46) 
(47) 
(48) 
(49) 



where E = [K 2 + q 2 sldc /A + m 2 } 1 ^ 2 . As in the previous case of the fireball geometry (Appendix 
A), one has to calculate X\ and X 2 - While the latter vanishes again by symmetry, the former 
reads after using eq. (3.547.4) of in the ^-integration 



47T / dC ^(0^(01-^^1 

+KK 



E COsh T) r 

' E cosh r] r 
f 



h(a,b) 



(50) 



where a = K sinh^/T, b = q S ide r f, and the functions Iq, I\ (which have to be evaluated 
numerically) are defined by 



J (0,6) 
h(a,b) 



7T JO 
1 



7T JO 



cosh [a cos (j)} cos [b sin 0] , 

cos (ft sinh [a cos <ft] cos [b sin 0] = 



di (a,b) 
da 



(51) 
(52) 



It is obvious (see eqs. (8.411.1, 8.431.4) in pESfl ) that J is related to the Bessel functions 

JoM) = J (o) , J o (0,6) = J (6) . (53) 
The final result for the side-correlation function reads (as in the preceding Appendix) 



C*2,side — 1 ~\-X 2 /Xq , 

with Xi from (^D|) and X from fl4"5j ) with the single inclusive spectrum 
For the out-correlation the choice of momenta is 



(54) 



Pi 

u 
P2 



(E u K + q out /2,0,0) , 
(E 2 , K - g out /2, 0, 0) , 
(K°,K,0,0) , 
(Si - E 2 , g ut, 0, 0) , 



(55) 
(56) 
(57) 
(58) 
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where E\$ = [(K ± g ou t/2) 2 + m 2 ] 1 / 2 , K° = (E\ + E 2 )/2. As in the previous cases, the 77- 
and 0-integrations separate, and the final result for X\ and Z 2 reads 



Xi = 4^^ dC ry(C) 7-/(0 [/Cx(a,/5) JoM) + /?) Jo (a, 6) 

[/C (a, /?) Ji(a, 6) + /C («, /3) Ji(a, 6) 
J 2 = 47T^ 1 dCr / (C)r / (C)|-K° [t x {a^) J (a,b) - ^(a, (3) J (a,b) 

+K [/C (a, /3) Ji(a, 6) - /C (a, /3) Ji(a, 6) 



dC 
dry) 

dC J 

dry 

dC 

dr/1 
dCj 



(59) 



(60) 



where a = K° cosh r/ r /T, /3 = — E^Tf, a = K sinh r] r /T, b = q out rf, and (see eqs. 
(8.431.4, 8.431.5, 8.432.1, 8.476.4) in 



Ko(a,P) 
/Ci(a,/3) 

£i(a,P) 
Jo{a, b) 
Jx{a,b) 
Jo{a, b) 



dr] cos [/3 coshr/] e~ a cosh " = Re K (a - i(3) , (61) 

d V coshr; cosLS coshr/] e~ a coshr? = Re^i(a - i(3) = f 62 ) 

Jo aa 

/ dr/ sin [3 cosh 77] e~ a cosh " = Im K (a - i/3) , (63) 
jo 

f°dr/ coshr/ sin[/3 coshr/] e - QCOsh ' 7 = ImrY^a-^) = ? ( 64 ) 

Jo aa 

d0 cos [6cos0] cosh [a cos0] = Re Jo(a + ib) , (65) 

d0 cos0 cos [6cos0] sinh [a cos</>] = Re/i(a + ib) = ^^ Q r [ Q ' - , (66) 



71 JO 
1 



7T JO 
1 



da 



7T Jo 
1 



sin [b cos 0] sinh [a cos 0] = Im J (a + z6) 



(67) 



7T JO 



d(j) cos0 sin [6cos0] cosh [a cos</>] = Im/i(a + ib) = ' ^ • (6£ 



5a 



The final result for the out-correlation reads 



a 



2, out 



1 + 



I?+J 2 



2b(pi,±) 2b(p2,±) 



(69) 



as before, cf. (H), but with X 1; X 2 from (|59"D and (|60f) and X from (|45|) with the single 
inclusive spectrum 
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Figure Captions: 



Fig. 1: (a) the entropy density divided by T 3 (in units of s c /T 3 ), (b) the energy den- 
sity divided by T 4 (in units of T c s c /T 4 ) as functions of temperature (in units of T c ), (c) the 
pressure (in units of T c s c ), (d) the square of the velocity of sound as functions of energy 
density (in units of T c s c ). The solid lines correspond to AT = 0, the dotted curves to 
AT = 0.1 T c . Quantities for the ideal gas equation of state (with dn degrees of freedom) 
are represented by dashed lines. The ratio of degrees of freedom in the QGP to those in the 
hadronic phase is dq/dn = 37/3. The critical enthalpy density is T c s c ~ 0.75 GeV fm -3 for 
the case g?q = 37, du = 3. 

Fig. 2: The expansion of (a,b) a fireball and (c,d) a Bjorken cylinder (with initial time 
t = 0.1 R) for an ideal gas equation of state p = e/3. (a,c) temperature (in units of the 
initial temperature T ), (b,d) laboratory energy density (in units of the initial pressure p ) as 
functions of radial distance r from the origin (in units of the initial radius R of the system). 
The profiles are for times t = r + nXR, n = 0, 1, 2, 5, A = 0.99. Solid lines are obtained 
via the semi-analytical method of characteristics ||28|| , dotted curves are produced with the 
relativistic HLLE algorithm modified with Sod's operator splitting method. 

Fig. 3: (a,c,e) temperature profiles for the fireball expansion for times t = OAnXR, n = 
0, 1, ...,8 and an initial energy density eo = §T c s c / {dq/du + 1) = 0.1125 T c s c = e#. The 
profiles are alternatingly shown as full and dotted lines in order to better distinguish them. 
(b,d,f) show isotherms in the corresponding space-time diagrams. The isotherms are labelled 
with the corresponding temperatures in units of T c . Figs. (a,b) are calculated with AT = 0, 
(c,d) with AT = 0.1 T c , and (e,f) with the ideal gas equation of state. 

Fig. 4: Same as in Fig. 3, for e = 1.875 T c s c ~ eq. Profiles in (a) are for times t = 2n\R, 
in (c) for t = nXR, and in (e) for t = 0.5 nXR, n — 0, 1, 10. 

Fig. 5: Same as in Fig. 3, for eo = 18.75 T c s c ~ 10 €q. Profiles in (a,c,e) are for times 
£ = 0.6nA#, n = 0,l,...,10. 

Fig. 6: Purely longitudinal expansion in the Bjorken model for an initial energy density 
eo = 10T c s c . (a) energy density (in units of T c s c ), (b) entropy density (in units of s c ), (c) 
pressure (in units of T c s c ), (d) temperature (in units of T c ) as functions of (proper) time r 
(in units of r ). Solid lines are for AT = 0, dotted for AT = 0.1 T c , dashed for the ideal gas. 
Different cooling laws r a are indicated. 

Fig. 7: Same as in Fig. 3, but for the Bjorken cylinder expansion with To = 0.1 R. Profiles 
in (a,c,e) are for times t = tq + 0.1 nXR, n — 0, 1, 8. 

Fig. 8: Same as in Fig. 7, for an initial energy density e = 1.875 T c s c ~ €q. Profiles 
in (a,c,e) are for times t = r + 0.3 nXR, n = 0,1, 10. 
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Fig. 9: Same as in Fig. 7, for an initial energy density eo = 18.75 T c s c ~ 10 eg. Pro- 
files in (a,c,e) are for times t = r + 0.6 nXR, n = 0, 1, 10. 

Fig. 10: Lifetimes (in units of the initial radius R of the system) in the one-dimensional 
(Landau) expansion (cf. also |14|]) as a function of the initial energy density (in units of 
T c s c ). (a,b) are for cLq/cIh = 37/3, (c,d) for dg/dn = 3. The thick lines in (a,c) are for 
AT = 0, in (b,d) for AT = 0.1 T c . Thin lines correspond to the ideal gas case. Solid lines 
are for T = 0.7 T c , dotted for AT = 0.9 T c , dashed for T = T c . 

Fig. 11: The same as in Fig. 10, but for the fireball expansion. 

Fig. 12: The same as in Fig. 10, but for the Bjorken cylinder expansion with r = 0.1 R. 
Fig. 13: The same as in Fig. 12, but for r = 1/3 T . 

Fig. 14: Two-pion correlation functions in (a,c,e) side- and (b,d,f) out-direction for the 
fireball expansion with eo = 1.875 T c s c . The average pion momentum is K = 300 MeV, the 
initial radius R was fixed to be 5 fm. The critical temperature was taken as T c = 160 MeV, 
the pion mass is m = 138 MeV. (a,b) are for AT = 0, (c,d) for AT = 0.1 T c , (e,f) for the 
ideal gas. Solid, dotted, and dashed lines are the correlation functions as calculated along 
the T = 0.7 T c , 0.9 T c , and T c isotherm, respectively. 

Fig. 15: The same as in Fig. 14, but for the Bjorken cylinder expansion with r = 0.1 R and 
eo = 18.75 T c s c . The average transverse pion momentum is K = 300 MeV, the longitudinal 
momenta of the pions vanish. 

Fig. 16: The same as in Fig. 11, but for the ratio i? out /-Rside- 
Fig. 17: The same as in Fig. 12, but for the ratio -R out /-Rside- 
Fig. 18: The same as in Fig. 13, but for the ratio -R ou t/-Rside- 
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